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1. Introduction 

The aim of this note is two-fold. First, we prove a conjecture of [5, 2] concerning 
the existence of a recurrence in k > satisfied by the integrals 

1 f°° f°° dx\dx<i ■ ■ ■ dx n 



o Jo 



(coshxi + • • • + coshx„) fe+1 ' 



O . for n = 1,2,.... These integrals naturally arose during the analysis of parts of the 

Ising theory of solid-state physics [5] . In [2] only the first four cases of Theorem 1 
f\J . below were proven and the proofs relied on the ability to express the corresponding 

^ ' integrals in (3) below as Meijer G-functions, something which fails for n > 4. 

A second aim is to advertise the power of current symbolic computational tools 
and related algorithmic developments to settle such questions. For this reason we 
give quite detailed exposition of the methods entailed. 

Our main result (Theorem 1) is better phrased in terms of 

O; _n!T(fc + l) 

C„.fe .— — ■ ^n,k- IT; 

o" 

Theorem 1 (Linear Recursion). For any fixed n 6 N, the sequence c Uy k enjoys a 
linear recurrence with polynomial coefficients of the form 

(k+l) n+1 C n , k + J2 Pn,j( k ) C n,k +j =0, (2) 

2<j<n 
j even 

with degP n j < n + 1 — j. 

Substituting (1) and simplifying by n\T(k + l)/2" yields 

(k + i) n+1 c n , k + ]T p nj (k)(k+j)(k + j-i)---(k + i)a hk+3 =o, 

2<j<n 
j even 

which is [2, Conjecture 1], with extra information added on the origin of the linear 
factors for the recurrence in C n k ■ 
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The starting point of our proof is the integral representation [2, Eq. (8)] 

c n .k = / t k K (t) n dt, (3) 
Jo 

where K is the modified Bessel function on which much information is to be found 
in [1, Ch. 9]. The key properties of Kq that we use are as follows. First 

/>OC 

K {t) = / e- tcosh ^dx 
Jo 

which explains how the integrals in (3) arise. Moreover, we have 

- a linear differential equation: (8 2 — t 2 )Ko(t) = with 8 := td/dt; 

- the behaviour at the origin: Ko(t) ~ — lnt, t — * 0; 

- and behaviour at infinity: Ko(t) ~ \J n/2te~ t , t — > +oo. 

These last two properties show that the integral (3) converges for any complex k 
subject to 5Jfc > —1 where it defines an analytic function of k. The recurrence 
of Theorem 1 then gives the integral a meromorphic continuation to the whole 
complex plane with poles at the negative integers. 

2. Existence of a Recurrence 

The theory of D-finite functions leads to a direct proof of the existence of a 
recurrence such as (2) in a very general setting, together with an algorithm. 

Recall that a function is called D-finite when it satisfies a linear differential 
equation with polynomial coefficients. A good introduction to the basic properties 
of these functions is given in [14]. What makes these functions appealing from the 
algorithmic point of view is that they live in finitc-dimcnsional vector spaces and 
thus many of their properties can be explicitly computed by linear algebra in finite 
dimensions. In this setting, the following proposition is easily obtained. It is a 
generalization of our main Theorem 1, except for the absence of degree bounds. 

Proposition 1. Assume that f(t) obeys a homogeneous linear differential equation 

a r (t)f^(t) + --- + a o (t)f(t)=0, 

with polynomial coefficients ai(t) in C[t\. For a fixed n G N\ {0}, let T be a path 
in C such that for any k,j £ N the integrals 

h,j ~ J^ k {f{t) n ) ij) dt 

converge and the limits of the integrand at both endpoints coincide. Then the se- 
quence {Ik,o} obeys a linear recurrence with coefficients that are polynomial in n 
and k and which can be computed given the coefficients 's. 

Wc give the proof in two steps. The first one is classical and can be found for 
instance in [14, Thm.6.4.9]. 

Lemma 1. D-finite functions form an algebra over the rational functions. 

This means that any polynomial in D-finite functions with rational function 
coefficients defines a functions that is itself D-finite. In particular Kq satisfies a 
linear differential equation. 
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Proof. The proof is effective. The difficult part is the product. The derivatives 
of two D-finitc functions / and g live in finite-dimensional vector spaces gener- 
ated by /,..., /' r_1 ) and g, . . . , g^ -1 - 1 . Therefore by repeated differentiation the 
derivatives of a product h := fg can be rewritten as linear combinations of the 
terms f {i) g u \ < i < r, < j < s which generate a vector space of dimension 
at most rs. It follows that the rs + 1 successive hS k \ k = 0, . . . ,rs are linearly 
dependent. A linear dependency between them can be found as the kernel of the 
linear map (Ao, . . . , A rs ) i— ► XqIi + • • • + X rs h^ rs \ Any such linear dependency is a 
linear differential operator annihilating fg. □ 

The corresponding algorithm is implemented, among other places, in the Maple 
package gf un [13]. 

Example 1. Here is how the function gf un [poltodif f eq] is invoked to compute 
a differential equation for Kq. 

> eqK0:=t*diff (t*diff (y(t) ,t) ,t)-t~2*y(t) ; 

d 2 d 
eqKO := t{t—y{t) + -y{t)) - t 2 y(t) 

> gfun [poltodif feq] (y(t)"4, [eqKO] , [y(t)] ,y(t))=0; 

t^yit) + 10< 3 ^yW - (20t 4 - 25t 2 )^y(0 

d 2 d 
- (120t 3 - 15t)—y(t) + (6At 4 - 152t 2 + 1)— y(t) + (128t 3 - 32t)y(i) = 
dt z dt 

Example 2. Here are the corresponding steps of the calculation for the smaller 
example Kq-. 

h = Kl, ti = 2K K , h" = 2K'q 2 - 2t- 1 K K l Q + 2K$, 

h (3) = -Qt- l K 12 + 4(2 + r 2 )K K' - 2t- 1 K 2 , 

where whenever possible we have replaced Kq by Kq — t~ 1 K{ ) . Then we find the 
vector (— 4i, 1 - At 2 , 3t, t 2 ) in the kernel of 

'10 2 -2* -1 \ 

2 -It' 1 4(2 + t- 2 ) . 
K 2 -6*- 1 / 

This vector then produces a differential equation satisfied by Kq: 

t 2 y^ + 3ty" + (1 - At 2 )y' - Aty = 0. 

Proof, (continued) The second step of the proof of Proposition 1 starts by expand- 
ing the differential equation for h := /" as 



^ </, ,/'// = 0, 



for scalars dij. This is then multiplied by t k and integrated along T. Use of the 
convergence hypotheses then allows us to deduce that 



t k+l h^ dt = 0. (4) 
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Now, integration by parts gives 

t k+l h {j) dt = t^h^-V 



-(k + i) / t^' 1 ^'^ dt, 



o 

= (-k- i){-k -i+ 1) ■■■(-k-i+j - l)h+i-j, 

the last equality following by induction. Adding the contributions of all the terms 
in (4) finally yields the desired recurrence over Ik- □ 

Example 3. For h := Kq, the computation gives 

/ t k+2 h^ + 3t k+1 h" + [t k - 4t k+2 )ti - 4t k+1 hdt = 0, 
Jo 

whence 

{-k - 2)(-k - l)(-k)c 2 , k -l + 3(-fc - l)(-*)C2,fc-l 

+ (-k)c 2>k _! - 4(-fc - 2)c 2 . k+1 - 4c 2:fe+ i = 0. 

Once simplified, this reduces to 

4{k + l)c 2 , k +i =fc 3 c 2 ,fc-i. (5) 

Example 4. Quantum field theorist David Broadhurst has recently studied [6] the 
vacuum- diagram integrals for n > 0: 

/•oo 

V(n,a,b):= / x 2n+1 K%{x){xK' n (x)) b dx 
Jo 

and provides the recursion 

2(n+l)V(n,a,b) + aV(n,a-l,b+l) + bV(n + l,a+l,b-l) = 

which preserves N := a + b > and allows one to reduce to V values with ah = 0. 
Note that K = —K±. Proposition f applies to n i— > V(n, a, b) for each a and b and 
as in Section 4 below leads to very efficient code for the recursion. The difficult 
question of understanding the initial values is discussed in [6] and [3] . 

Mellin transform. As the proof indicates, Proposition f is not restricted to integer 
values of k. In particular, the method gives a difference equation for the Mellin 
transform 

r + co 

h*(s) := / t^h^dt 
Jo 

provided the appropriate convergence properties are satisfied. This difference equa- 
tion then gives a meromorphic continuation in the whole complex plane. The most 
basic example is T(s): starting from the elementary differential equation y' + y = 
for h{t) = exp(— t) leads to the classical functional equation T(s + 1) = sT(s). 
Coefficients. The path T can also be a closed contour. For instance, if h is analytic 
at the origin, then the fcth coefficient of its Taylor series at the origin is given by 
the Cauchy integral 

— / 4tt dt, 

2iri J t k+1 ' 

where the contour encloses the origin and no other singularity of h. The algebraic 
manipulations are the same as in the previous case, followed by replacing k by — fc — f 
and the sequence c k by the sequence U-k-i- 
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For instance, if we apply this transform to the functional equation for T, we get 
first c_s = —sc-s-i and then u s -i = —su s , which is the expected recurrence for 
the sequence of coefficients u s = (— l) s /s! of exp(— t). 

Similarly, starting from (5) we obtain the mirror recurrence 

4fc5 fc _i = (k+l) s c k +i- 
Observe that this is obeyed by the coefficients of In 2 (t) in the series expansion 

Kg® - ,»>(«> (1 + 1,* + ^ + ^ + ...) 

+ ln(t) ^2 7 -21n2 + ( 7 - ln2-i)t 2 + ...^) + ((ln2- 7 ) 2 + ...), t- 

The Frobenius computation of expansions of solutions of linear differential equations 
at regular singular points (see, e.g., [11]) explains why this is so. 
Further generalizations. The sequence {t k } in the integral of the proposition 
can be replaced by more general sequences of functions that satisfy both a linear 
recurrence in k and a linear differential equation in t. Provided that proper analytic 
conditions are satisfied at the endpoints of the path, the same result will hold. This 
makes it possible to compute, for instance, recurrences for the Fourier coefficients 
with respect to various bases such as orthogonal polynomials, Bessel functions,. . . . 

The algorithm that applies in this case is called creative telescoping as discovered 
by Zeilberger [16] and further automated in [9, 8]. Again, the computation boils 
down to linear algebra in a suitably constructed finite-dimensional vector space. 

In summary, all these algorithms succeed in making effective and efficient the 
familiar method of differentiation under the integral sign and integration by parts. 

3. Proof of the Main Result 

If A is a linear differential operator, the operator of minimal-order annihilating 
the nth power of every solution of A is called its nth symmetric power. Because 
of its role in algorithms for differential Galois theory [15] there has been interest 
in efficient algorithms computing symmetric powers. In the case of second order 
operators, such an algorithm has been found in [7]. We state it in terms of the 
derivation 8 := td/dt in order to get better control over the coefficients of the 
resulting recurrence — but the statement and proof hold for any derivation. 

Lemma 2 (Linear Differential Equation, [7]). Let A = 8 2 +a(t)9 + b(t) be a linear 
differential operator with rational function coefficients a(t) and b(t). Let Lq = 1, 
L\ = 6 and for k = 1, 2, . . . , n define the operator Lk+i by 

L k+1 := (0 + ka)L k + bk(n -k+ l)L k _ v (6) 
Then, for k = 0, . . . , n + 1, and for an arbitrary solution y of Ay = 0, 

L k y n = n(n - 1) • ■ • (n - k + l)y n - k {6y) k 

and in particular L n+ iy n = 0. 

[This recursion can be viewed as an efficient computation of the kernel which 
was described in the previous section, taking advantage of the special structure of 
the current matrix.] 
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Proof. The proof is a direct verification by induction. For k = and k = 1 the 
identity reduces respectively to y n = y n and 6y n = ny n ~ 1 8(y) which are obviously 
true for any function y. Assuming the identity to hold up to k > 1, the heart of 
the induction is the rule for differentiation of a product 9(uv) = 6(u)v + u9(v): 

e{y n ' k {6y) k ) = e(y n - k ){6y) k +y n - k 6{{ey) k ) 

= (n - k)y n - k -\6y) k+1 + ky n - k {e y ) k -\6 2 y) 

= (n - k)y n - k - 1 {9y) k+1 + ky n - k {6y) k - 1 {aBy + 6y). 

Reorganizing terms concludes the induction. □ 

Example 5. In the case of Kg, we have a = and 6 = — t 2 . For n = 4, starting 
with Lo = 1 and L\ = 9, the recurrence of Lemma 2 gives 

L 2 = 6 2 - U 2 , 

L 3 = 6 3 -10t 2 9-8t 2 , 

Li = 6 4 - 16t 2 9 2 - 28t 2 9 + 8t 2 (M 2 - 2), 

l 5 = e 5 - 20t 2 e 3 - 60t 2 e 2 + 8t 2 (8t 2 - 9)6 + m 2 {At 2 - 1). 

The operator L5 annihilates Kq. It is a rewriting in terms of 9 of the equation of 
Example 1. 

Some of the patterns that emerge on this example can be proved in the general 
case. 

Lemma 3 (Closed Form). With the same notation as in Lemma 2, when A = 
2 — t 2 , Lfc may be written as 

fc-2 

L k =6 k +Y j af\t)e^ 

3=0 

where each is a polynomial in t 2 , divisible by t 2 and dega^ < k — j. 

Proof. Again, the proof is by induction. For k = and k = 1 we recover the 
definition of Lq and L\. For larger k, the recurrence (6) simplifies to 

L k+1 := 6L k - k(n -k+ l)t 2 L k ^. 

If the property holds up to k > 1, then this shows that the degree of L k +i in 8 
is k + 1, with leading coefficient 1 and also that the coefficient of 8 k in L k +i is 0. 
Extracting the coefficient of Qi then gives 

(fc+i) = faf-i + ( a f ') ^ k ( n - k + !)* 2 4* < j < fc - 2, 
° J \-k(n-k+l)t 2 , j = k-l. 

These last two identities give the desired degree bound and divisibility property for 
the coefficients a^ 1 ', < j < k — 1. □ 

We may now complete the proof of the main result. 

Proof, (of Theorem 1.) Lemma 3 shows that the coefficients of L n+ i can be rewrit- 
ten as 

L n+1 = 6 n+1 + *QMi (7) 

2<j<n 
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where the polynomials Qj satisfy deg Qj < n + 1 — j. 

Thanks to the properties of Kq recalled in the Introduction, an integration by 
parts yields 

r+oo 

/ t k+ i6 m {K%{t)) dt=(-l-k- j) m c n . k+J , (8) 
Jo 

for each m. Now, we multiply L„ + i/Tq from (7) by t k and integrate from to 
infinity: 




j even 



Integrating term by term and using (8) finally gives the recurrence 

(-k-l) n+1 c n , k + Qi{-l-k-3)Cn,k+i =0. 

2<j<n 
j even 

which is the desired one, up to renaming and sign changes. □ 

4. Algorithm 

In summary, we have a relatively straightforward algorithm to compute the linear 
recurrences for the c n ,k or C njk for given n. First, the operators L k can be computed 
as commutative polynomials L k as follows: 

dL 

L k+1 :=t—± + 9L k -k(n-k + l)t 2 L k _ 1 , l<k<n, (9) 
at 

with initial values Lq := 1 and L\ := 9. These polynomials L k coincide with the 
operators L k when the powers of 9 are written on the right of the monomials in t 
and 9. 

By collecting coefficients of t in we recover (7). Substituting — 1 — k — j 

for 9 in the coefficient of t 3 then produces the desired recurrence for c n _ k , while 
replacing c n _ k+ j by (k + 1) • ■ ■ (k + j)C n . k +j for all j produces one for C n . k . 

Example 6. We illustrate the process for n = 4. The last operator in Example 5 
may be rewritten as 

L 5 =9 5 - 4t 2 (56» 3 + 159 2 + 189 + 8) + 64£ 4 (6> + 2) 

and annihilates K^(t). Substituting — 1 — k — j for 9 in the coefficient of t 3 for 
j = 0, 2, 4 gives 

-(k + l) 5 c 4 , fc + 4(> + 2)(5fc 2 + 20k + 23)c 4 , fe+ 2 - 64(fc + 3)c 4 , fc+ 4 = 0. 
Since CA,h = §T(fc + l)Ci. k , this is equivalent to 

-|(fc+l) 4 C4,fc + 6(fc + 2) 2 (5fc 2 + 20fc + 23)C4,fc + 2-96(fc + 4)(fc + 3) 2 (fc + 2)C4,fe + 4 = 

which was proven by different methods in [2]. 

Here is the corresponding Maple code: 
compute_Q :=proc(n,theta,t) 
local k, L; 

L[0] :=1; L[l] :=theta; 
for k to n do 

L[k+1] : =expand(series ( 
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t*diff (L[k] ,t)+L[k] *theta-k* (n-k+1) *t~2*L [k-1] , 
theta, infinity) ) 

od; 

series (convert (L [n+1] ,polynom) ,t, infinity) 

end: 

rec_c : =proc (c : :name,n: :posint,k: :name) 
local Q, theta, t,j; 

Q:=compute_Q(n, theta, t) ; 

add(factor(subs(theta=-l-k-j ,coeff (Q,t, j)))*c(n,k+j) , j=0. .n+l)=0 

end: 

rec_C : =proc (C : :name,n: :posint,k: :name) 
local Q, theta, t, j , ell; 

Q:=compute_Q(n, theta, t) ; 
(-l)~(n+l)*(k+l)~n*C(n,k)+ 

add(f actor (subs (theta=-l-k-j , coef f (Q , t , j ) ) 

*mul(k+l+ell,ell=l. . j-1) ) *C(n,k+j ) , j=l . .n+l)=0 

end: 

On a reasonably recent personal computer, all recurrences for n up to 100 can be 
obtained in less than 5 minutes (further time could be saved by not factoring the 
coefficients). For example, the recursions for c 4i fc and C 4 .fc may be determined thus: 

> rec_c(c, 4, k) ; 

-(fe + l) 5 c 4 , fc + 4(fc + 2)(5fc 2 + 20fc + 23)c 4 , fe+2 - (64fc + 192)c 4 , fc+4 = 

> rec_C(C, 4, k) ; 

- (k + l) 4 C 4 , fe + 4(fc + 2) 2 (5fc 2 + 20k + 23)C 4 , fc+2 

- 64(k + 4)(k + 3) 2 (fc + 2)C 4 , fc+4 = 

The first six cases for C n ^ are 

= (fc+l)Ci I fc-(fc + 2)<7i, fc+2 (10) 

= (k + l) 2 C 2 . k -A{k + 2) 2 C 2 M + 2 (11) 

= (fc + l) 3 C 3 , fc -2(fc + 2)(5(fc + 2) 2 + l)C 3 , fc+2 

+9(fc + 2)(fc + 3)(A; + 4)C3 )fe+ 4 (12) 

= (fc+l) 4 C 4 , fc -4(A: + 2) 2 (5(A : + 2) 2 + 3)C 4 , fe+2 

+64(/c + 2)(/c + 3) 2 (/c + 4)C 4ifc+4 (13) 

= (/c + l) 5 C5,fc-(/e + 2)(35fc 4 + 280fc 3 + 882/s 2 + 1288fc + 73l)C 5:A; +2 
+ (k + 2)(fc + 3)(fe + 4) (259fc 2 + 1554k + 2435) C 5 ,fe +4 
-225(fc + 2)(fc + 3)(fc + 4)(k + 5)(fc + 6)C 5 , fe+6 (14) 

= (fc + l) 6 C 6 , fc -8(/c + 2) 2 (7A: 4 + 56A: 3 + 182fc 2 + 280A:+17l)C* 6 , fc+ 2 

+16(fc + 2){k + 3) 2 (k + 4) (49k 2 + 294fc + 500) C 6 , fc+4 

-2304(/c + 2)(k + 3)(k + 4) 2 (k + 5)(k + 6)C 6 , k +6- (15) 

as given in [2], but in which only the first four were proven (see also [12] for an 
earlier proof up to n = 4). Many more recursions were determined empirically 
using Integer Relation Methods — this relied on being able to compute the integrals 
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in (3) to very high precision — and led to the now-proven conjecture. The versions of 
these recurrences in terms of c n ^ instead of C n ,k were also determined empirically 
for n = 1, . . . , 6 in [10, Eqs. (lla-e)] for the enumeration of staircase polygons. 

Implicit in this algorithm is an explicit recursion for the polynomial coefficients of 
each recursion. In the case of (12) and (13) these recursions lead to new continued 
fractions for L_ 3 (2) and £(3) respectively [5, 2). These rely additionally on the 
facts that C 3 ,i = L_ 3 (2),C 3j3 = 2L_ 3 (2)/9 - 4/27 and C 4 ,i = 7C(3)/12, C 4>3 = 
7£(3)/288 — 1/48, [2]. Corresponding continued fractions arising from C 3i 2/C 3j o 
and 64,2/64,0 are determined in [3]. 



5. Another Example 

In [4] (to which we refer for motivation and references), the following "box inte- 
grals" have been considered 

B n (s)= f ... f {r 2 1 + ---+rly/ 2 dr 1 ---dr n , 
Jo Jo 

A„(s) = / . . . / ((n - qif H h (r„ - q n ) 2 ) s/2 dri ■ ■ ■ dr n dqi ■ ■ ■ dq n . 

Jo Jo 

As in the case of the C n> k we have considered here, a good starting point is provided 
by alternative integral representations for -Rs > 0: 



0Ferf(u) 



B n{~s)^—— \ U S -%u) n du, b{u) 

T{s/2) J a 2u 
A n(-s) = „, /0 , / u * l d{u) n du, d(u) = j — 



T( s /2) 

The first one is given explicitly as [4, (33)] and the second one can be derived sim- 
ilarly. From classical properties of the error functions, the functions b(u) and d(u) 
satisfy the linear differential equations 

ub"(u) + 2(1 + u 2 )b'{u) + 2ub(u) = 0, 
2u 2 d"'{u) + 4u(3 + u 2 )d"{u) + 4(3 + Au 2 )d'{u) + 8ud(u) = 0. 

This is exactly the set-up of our Proposition 1. We thus deduce the existence of 
linear difference equations (wrt s) for both B n and A n . The fast computation of 
the difference equation for B n follows directly from the Algorithm of the previous 
section, and for instance, we get 

(s + 9)(s + 10)(s + ll)(s + 12)B 4 (s + 8) - 10(s + 8) 2 (s + 9)(s + 10)B 4 (s + 6) 
+ (s + 6)(s + 8)(35s 2 + 500s + 1792)B 4 (s + 4) - 2(25s + 148)(s + 4)(s + 6) 2 B 4 (s + 2) 

+ 24(s + 2){s + 4) 2 (s + 6)B 4 {s) = 0. 

The recurrence holds for all s by meromorphic continuation. A result on the shape 
of this recurrence for arbitrary n could be obtained along the lines of Lemma 3. 

Acknowledgements The authors wish to express their thanks to David Broad- 
hurst for directing them to several relevant references and for his many incisive 
comments. 
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